Microscopic structure of liquid hydrogen: a neutron diffraction experiment 
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We have measured the center-of-mass structure factor S(k) of liquid para-hydrogen by neutron 
diffraction, using the D4C diffractometer at the Institute Laue Langevin, Grenoble, France. Fhe 
present determination is at variance with previous results obtained from inelastic neutron scattering 
data, but agrees with path integral Monte Carlo simulations. 



PACS numbers: 67.20.+k, 61.20.Ja, 61.12.-q 

Among the simple liquids, the microscopic structure 
factor, S(k), of liquid hydrogen has been one of the most 
jealously kept secrets of nature. This is a consequence of 
several experimental problems. Due to the small num- 
ber of electrons, an X-ray determination of the hydrogen 
structure factor is not an easy task. On the other hand, a 
neutron diffraction experiment is particularly difficult be- 
cause of the relevance of inelastic scattering events: their 
size, mainly determined by the ratio between the neutron 
and the nuclear mass jjj , makes the standard correction 
techniques unsuitable to hydrogen. Due to the doubled 
nuclear mass, the experiment analysis is less demand- 
ing for deuterium. Nonetheless, the first reliable neutron 
diffraction measurement of the microscopic structure of 
liquid deuterium is relatively recent and was carried out 
using a small-angle, time-of- flight, diffractometer (SAN- 
DALS, at ISIS, UK) that was specifically built for light 
mass liquids. Using this instrument, the structure factor 
of liquid D2 was measured in the vicinity of the triple 
point Q and close to the freezing transition ||. By 
combining the data of Ref. || with those measured in 
a further diffraction experiment, performed on a reac- 
tor source (7C2 diffractometer at Laboratoire Leon Bril- 
louin, Saclay, France) Q), an improved determination of 
the structure factor of liquid deuterium was obtained. 
This turned out to be in excellent agreement with path 
integral Monte Carlo (PIMC) simulation results S. 

The knowledge of S(k) for deuterium is not easily 
transferred to hydrogen, due to the different role played 
by quantum effects ||. On the other hand, for the case 
of hydrogen, the overwhelming ratio between the inco- 



herent and coherent neutron scattering cross section of 
the proton makes it extremely difficult to extract the in- 
termolecular response, which carries the structural infor- 
mation, from the large intramolecular contribution. The 
first successful attempt to obtain structural information 
on liquid hydrogen, carried out on SANDALS, allowed us 
to determine the thermodynamic derivatives of S(k) [Q. 
These were found in a rather good quantitative agree- 
ment with the PIMC simulation results || . 

Due to the well recognized experimental difficulties, al- 
ternative methods have been suggested to determine S(k) 
for liquid hydrogen. Bermejo and co-workers ^ have re- 
ported a determination of the structure factor of liquid 
para-hydrogen using the results of an inelastic neutron 
scattering experiment and the sum rule which relates the 
dynamic structure factor, S(k,u>), to the static one, S(k), 
PI: 



S(k) = 



duj S{k, uj). 



(i) 



The results were qualitatively reasonable, but not fully 
convincing on a quantitative basis p| . In fact, the main 
peak of S(k) appeared rather high, exceeding the value 
of 2.85 which marks, according to the Hansen- Verlet cri- 
terion fl^| , the onset of the freezing transition. Even 
though Hansen and Verlet formulated this rule for simple 
model systems (classical monatomic particles interacting 
through a Lennard- Jones potential) it is hard to believe 
that for hydrogen, a genuine quantum system for which 
one expects an overall broadening and damping of the 
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structural features, S(k) should reach such a large value 
in the liquid phase. 

More recently, the same method, i.e. integrating the 
dynamic structure factor at constant k, was applied by 
Pratesi et al. [[l3| using the X-ray inelastic scattering 
data obtained on the ID16 beamline at the European 
Synchroton Radiation Facility (Grenoble, France). Their 
results for S(k), taken at T = 31.5 K and molecular num- 
ber density n =21.5 nm -3 , were not much extended in 
k, barely exceeding the main peak position of the struc- 
ture factor, and with rather large error bars. However, 
the peak height turned out considerably lower than the 
neutron data of Ref. ^| and consistent with their PIMC 
simulation results. The X-ray diffraction measurement, 
carried out in parallel E3|, showed much smaller error 
bars and still a good agreement with the PIMC simula- 
tion. However, in this case, the accessed fc-range did not 
even reach the peak position of S(k). 

In this context, we have carried out a new experiment 
aiming to determine directly, for the first time by neutron 
diffraction, the center-of-mass structure factor of liquid 
para-hydrogen and the results constitute the object of the 
present letter. The measurements were carried out using 
the D4C diffractometer of the Institut Laue Langevin 
(Grenoble, France) in its standard configuration with an 
incident wavelength Ao = 0.6933 A. Of the four thermo- 
dynamic states investigated, we report here the results at 
T = 17.1±0.1 K andp = 29.9±0.1 bar (n = 22.95±0.05 
nm~ 3 [0). The liquid hydrogen sample was condensed 
directly in a vanadium cylindrical container (6 mm in- 
ternal diameter, 0.2 mm wall thickness) imbedded in the 
low pressure buffer gas of an orange cryostat. At the 
bottom of the vanadium container, out of the neutron 
beam, we had inserted some powder of a paramagnetic 
catalyst made of Cr 2 03 on an AI2O3 substrate, in order 
to speed up the conversion from ortho- to para-hydrogen. 
The relative concentration of the two species was moni- 
tored looking at the low momentum transfer portion of 
the diffraction pattern Q] . Twenty-four hours of conver- 
sion time were allowed to stabilize the sample. Then, 
the stability of the low-fc portion of the diffraction pat- 
tern confirmed that the sample had reached the ther- 
modynamic equilibrium concentration, calculated to be 
99.96% rich in para-hydrogen. In order to check the over- 
all stability of the experimental setup, as well as of the 
sample, for each thermodynamic point we took several 
measurements in repeated runs. 

In order to properly subtract the background contribu- 
tion from the raw data, we have carried out an additional 
measurement, without changing the experimental setup 
and filling the sample container with a small amount of 
gaseous 3 He. The density of the gas was chosen such as to 
match the scattering power of the liquid para-hydrogen 
sample with the absorption power of 3 He. The details 
of the procedure can be found in Ref. 
value <r(Ao) = 43.4 barn, obtained from Ref. 
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FIG. 1: Diffraction pattern of liquid para-hydrogen at T = 
17.1 K and p = 29.9 bar. What in the figure appears as a bold 
line is composed, instead, by the ensemble of the experimental 
points. The lack of intensity in the low-fc region is peculiar 
of the para species. The decay at high momentum transfer 
is determined by the finite energy of the incident neutrons 
(see text). The tiny structure in the region k ~ 2 A -1 is the 
signature of the main peak of the intermolecular structure 
factor S(k). The dashed line is the calculated self molecular 
part. 



molecular hydrogen scattering cross section at the inci- 
dent neutron wavelength. In Fig. 1, we show the diffrac- 
tion pattern of liquid para-hydrogen after subtraction of 
background and container scattering, properly corrected 
for attenuation. The diffraction spectrum is dominated 
by the self molecular part, which is produced by the in- 
tramolecular structure. The distinct coherent contribu- 
tion, containing the required information on the inter- 
molecular structure factor, appears as a tiny undulation 
superimposed on the much larger intramolecular struc- 
ture in the region k ~ 2 A -1 . The difference between 
the experimental and calculated patterns (see Fig. 1) is 
mainly due to multiple scattering effects. 

The data treatment needed to extract the structure 
factor requires the calculation of the self part of the 
molecular double differential neutron cross section of 
para-hydrogen. This can be done using the Young and 
Koppel (YK) theory p7| . Here, the intermolecular inter- 
actions are totally neglected and the hydrogen molecules 
are modeled as a set of non-interacting particles possess- 
ing their relevant internal degrees of freedom. Thus, each 
molecule is considered separately and the intramolecu- 
lar roto-vibrational modes, and spin correlations, are ex- 
plicitly taken into account. The vibrational modes are 
considered harmonic and the rotations are free. This 
model applies well to our case because, at low temper- 
atures, only the ground vibrational state is populated 
and, using thermal neutrons, no vibrational transition is 
allowed. Moreover, as long as the system is in the liquid 
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phase and is not subject to high pressures, the anisotropic 
components of the intermolecular potential are negligible 
p8|, 19, M. The vibrational-rotational coupling can be 



accounted for, in an effective way, using the experimental 
values for the roto- vibrational levels J2l| . Within the ap- 
plicability limits of this model, a rigorous calculation of 
the double differential neutron cross section is possible. 

It is well known |||] that for hydrogen, due to quan- 
tum effects, the average center-of-mass kinetic energy is 
different from the classical value (E^) = (3/2)/cbT and 
depends on density fl23fl . This quantity can be directly 
measured from an inelastic neutron scattering experi- 
ment [Q. We have shown |l(| that a modified Young 
and Koppel (MYK) model can be defined to account 
for this property, using an effective temperature in the 
width of the ideal-gas Gaussian that describes the dy- 
namic structure factor of the molecular centers of mass. 
By means of the MYK model, the self portion of the 
double differential neutron cross section could be inte- 
grated, at each scattering angle 8, over the energies of 
the scattered neutrons. The result, properly normalized 
as described below, is represented by the dashed line in 
Fig. I. 

As far as the multiple scattering effects are concerned, 
we carried out a Monte Carlo (MC) simulation using, 
as input, the calculated MYK cross section. We found 
that the multiple scattering contribution to the cross sec- 
tion is a rather smooth function of k. However, since 
its intensity is larger than the intermolecular part (see 
Fig. 1), the accuracy of a typical multiple scattering cal- 
culation might be insufficient for a reliable extraction of 
S(k). On the other hand, the similarity between the mea- 
sured diffraction pattern and the calculated intramolec- 
ular cross section suggests to consider an overall, rather 
unstructured, generalized background, that includes the 
multiple scattering and should be represented by a low- 
order polynomial function P(k). Therefore, instead of at- 
tempting a difficult evaluation of the multiple scattering, 
we applied a different procedure where the order and co- 
efficients of the polynomial, as well as the overall normal- 
ization constant, were obtained from a fit. A model func- 
tion, including P(k) and the calculated intramolecular 
diffraction pattern, was fitted to those data points mea- 
sured at k— values where an independent knowledge of 
the intermolecular structure factor is available: namely, 
we used S(k) ~ 5(0) at low k, as is typical of a dense 
liquid (with S(0) related to the isothermal compressibil- 
ity), and S(k) ~ 1 for large k. The /c-range that contains 
the experimental information on the sought-for structure 
factor was excluded from the fit. We found that a fifth 
order polynomial is sufficient to give a satisfactory de- 
scription of the generalized background function which 
turned out to be in a semiquantitative agreement with 
the MC evaluation of the multiple scattering intensity. 

The fitted polynomial has been subtracted, together 
with the calculated intramolecular contribution, from the 
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FIG. 2: Measured intermolecular (i.e. center-of-mass) struc- 
ture factor of liquid para- hydrogen at T = 17.1 K and 
p = 29.9 bar. The error bars are determined as described 
in the text. Data are rebinned every five points in the region 
where S(k) is weakly k— dependent (i.e. for k < 1.5 A -1 and 
k > 2.8 At 1 ). The full line represents the results of a PIMC 
simulation using the Norman et al. ]27j spherical potential. 
The value for S(k = 0) is obtained from the isothermal com- 
pressibility pH |. 



experimental diffraction pattern. The result gives, after 
normalization, the intermolecular term in the measured 
cross section that can be written as: 



dn) 



= A SlSC (k)u(k)[S{k) - 1] 



(2) 



coll 



where A SjSC is the Paalman and Pings [g5|] attenuation 
coefficient and u(k) is the molecular form factor f26| . It 
is worthwhile to note that, in the data analysis, we have 
taken into account self absorption as well as the finite 
dimensions of sample and detectors. The final result for 
S(k) is reported in Fig. 2. The error bars account for 
both the propagated statistical errors from the raw spec- 
trum and an estimate of the systematic errors coming 
from the fitting procedure. For example, we have taken 
into account, among other things, the effect of the vari- 
ation of the fitting intervals in the determination of the 
polynomial background. In the same figure, we com- 
pare the present experimental neutron diffraction results 
with the corresponding quantity obtained from a PIMC 
computer simulation JIl[ using the Norman, Watts, and 
Buck |27[] (NWB) intermolecular potential. The simula- 
tion conditions were T = 17.1 K and n = 23.0 nm -3 . 
The agreement is not perfect but still is rather good and 
both the experimental and simulation data consistently 
show no trace of an anomalously high main peak in S(k), 
as it might have been inferred from Ref. |9[ 

We compare the present results with a corresponding 
analysis carried out on liquid deuterium, in similar ther- 
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modynamic conditions and using the same intermolec- 
ular potential (NWB). In that case, the agreement be- 
tween simulation and experiment was extremely good, 
and data were of such a quality to allow a good dis- 
crimination among different intcrmolecular interaction 
functions, namely NWB and Lennard-Jones In the 
present case, an overall quantitative agreement with sim- 
ulation is also found, although the experiment is much 
more difficult. Consequently, the resulting uncertainties 
are larger, and do not allow us to clearly discriminate 
among various choices of the intermolecular potential. 

In order to be more specific, we examined the struc- 
tural information that is obtained from the simulations, 
using two different potential functions. To this aim, we 
compared the present simulation data at T = 17.1 K and 
n = 22.22 nm~ 3 , based on the NWB potential, with ex- 
isting PIMC data at T = 17.2 K and n = 22.1 rxnT 3 
p8| , obtained using the Silvera and Goldman (SG) in- 
termolecular potential . The only observed difference 
was a lower main peak, by ~ 4% in [S(k) — 1], for the 
simulation results based on the SG potential function. 
Considering that temperature and density differ, in the 
two simulations, by ~ 0.6% and ~ 0.5% respectively, we 
can mainly attribute to the different choice of the inter- 
action potential the found deviation of ~ 4%. This, in 
turn, is a little smaller than the present size of the error 
bars in the main peak region. 

Finally, on closer inspection of Fig. 2, data seem to 
suggest a slightly narrower main peak than in the simu- 
lation. We exclude that such an effect could be assigned 
to the particular choice of the intermolecular potential. 
In fact, by comparing the NWB and the SG simulations, 
no observable difference in the main-peak width could 
be evidenced. Of course, we could attribute this effect 
to some remaining systematic experimental error that 
we were unable to get rid of. However, it is interest- 
ing to note that a similar qualitative result, concerning 
the width of S(k) , was also found in Ref. ^ and attributed 
to the presence, in the real liquid, of significantly longer 
range correlations than those depicted by the simulation. 
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